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Abstract 

The results and conclusions by Ponomarev and Vdovin [Phys. Rev. C 72, 034309 (2005)] are 
inadequate to judge the applicability of the modified BCS because they were obtained either in the 
temperature region, where the use of zero-temperature single-particle spectra is no longer justified, 
or in too limited configuration spaces. 
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The modified BCS theory (MBCS) was proposed and developed in as a micro- 

scopic approach to take into account fluctuations of quasiparticle numbers, which the BCS 
theory neglects. The use of the MBCS in nuclei at flnite temperature T washes out the 
sharp superfluid-normal phase transition. This agrees with the predictions by the macro- 
scopic theory ^, the exact solutions |5], and experimental data j^. The authors of [3] 
claimed that the MBCS is thermodynamically inconsistent and its applicability is far below 
the temperature where the conventional BCS gap collapses. The present Comment points 
out the shortcomings of 0]. We concentrate only on the major issues without repeating mi- 
nor arguments already discussed in or inconsistent comparisons in Fig. 9 and footnote 
[11] of (See B). 

1) The application of the statistical formalism in flnite nuclei requires that T should be 
small compared to the major-shell spacings (~ 5 MeV for ^^°Sn). In this case zero-T single- 
particle energies can be extended to T 7^ 0. As a matter of fact, the T-dependent Hartree- 
Fock (HF) calculations for heavy nuclei in j^l have shown that already at T > 4 MeV the 
effect of T on single-particle energies cannot be neglected. We carried out a test calculation 
of the neutron pairing gap for ^^''Sn, where, to qualitatively mimic the compression of the 
single-particle spectrum at high T as in j^, the neutron energies are e'j = ej(l + 7T^) with 
7 = —1.2 X 10"^ if \j) < |l5'9/2)- For \j) above |l5'9/2), we took 7 equal to 0.49 x 10"^ 
and —0.7 x 10~^ for negative and positive ej, respectively. The obtained MBCS gap has a 
smooth and positive T dependence similar to the solid line in Fig. 7 of ^ | with a flat tail of 
around 0.2 MeV from T = 5^MeV up to T = 7 MeV. For the limited spectrum used in the 
calculations of Ni isotopes [2], the major-shell spacing between (28-50) and (50-82) shells 
is about 3.6 MeV, so the region of valid temperature is T <^ 3.6 MeV. Hence, the strange 
behaviors in the results obtained at large T for ^^°Sn and Ni isotopes in ^ occurred because 
the zero-T spectra were extended to too high T. Moreover, the conflguration spaces used 
for Ni isotopes are too small for the MBCS to be applied at large T. The same situation 
takes place within the picket-fence model (PFM) analyzed below. 

2) The virtue of the PFM is that it can be solved exactly in principle at T = 0. However, at 
T 7^ the exact solutions of a system with pure pairing do not represent a fully thermalized 
system. As a result, temperatures deflned in different ways do not agree Q|. The limitation 
of the conflguration space with = 10 causes a decrease of the heat capacity C at Tm > 1.2 
MeV(Schottky anomaly) ^ (See Fig. 4 (c) of 7]). Therefore, the region of T > 1.2 MeV, 
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FIG. 1: (a) MBCS quasiparticle- number fluctuations 6J\fj within the PFM versus single-particle 
energies at several T. Lines connect discrete values to guide the eyes; numbers at the lines show 
the values of T in MeV; (b) BCS and MBCS gaps for iV = 10 and 17 = 11 (G = 0.4 MeV). 

generally speaking, is thermodynamically unphysical. The most crucial point here, however, 
is that such limited space deteriorates the criterion of applicability of the MBCS (See Sec. 
IV. A. 1 of 3]), which in fact requires that the line shapes of the quasiparticle- number 
fluctuations 5M± = \/ nj{l — rij) should be included symmetrically related to the Fermi level 
[Fig. 1 (f) of pI is a good example]. The dashed lines in Fig. ^ (a) shows that, for = 
10 particles and n = 10 levels {G = 0.4 MeV), at T close to 1.78 MeV, where the MBCS 
breaks down, 5A/^- are strongly asymmetric and large even for lowest and highest levels. At 
the same time, by just adding one more valence level {Q = 11) and keeping the same A^ = 
10 particles, we found that SN'j are rather symmetric related to the Fermi level up to much 
higher T [solid lines in Fig. ^ (a)]. This restores the balance in the summation of partial 
gaps 6Aj [3]. As a result the obtained MBCS gap has no singularity at < T < 4 MeV [Fig. 
^(b)]. The total energy and heat capacity obtained within the MBCS also agree better with 
the exact results than those given by the BCS [Fig. E]. It is worth noticing that, even for 
such small A^, adding one valence level increases the excitation energy E* by only ~ 10% at 
T = 2 MeV, while at T < 2 MeV the values of E* for = 10 and 11 are very close to each 
other. We also carried out the calculations for larger particle numbers A^. This eventually 
increases Tm, and also makes the line shapes of SJVj very symmetric at much higher T. For 
= 50 and 100, e.g., we found Tm > 5 MeV, and the MBCS gap has qualitatively the same 
behavior as that of the solid line in Fig. ^ (b) up to T ~ 5 - 6 MeV. However, for large A^ 
the exact solutions of PFM turn out to be impractical as a testing tool for T ^ 0. Since all 
the exact eigenstates must be included in the partition function Z, and, since for A^ = 50 
e.g., the number of zero-seniority states alone already reaches lO^'', the calculation of exact 
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FIG. 2: Total energies (a) and heat capacities (b) within the PFM for (N = 10, 0, = 11, G = 
0.4 MeV) versus T. Dotted, thin-, and thick-sohd hnes denote the BCS, MBCS and exact results, 
respectively. A quantity equivalent to the self-energy term —GY^^vj, not included within BCS 
and MBCS, has been subtracted from the exact total energy. 
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FIG. 3: hj (a) and Cj (b), obtained within BCS for 5 lowest levels in the PFM with = 10 versus 
T. In (a) the solid and dashed lines represent hj and quasiparticle energies Ej, respectively. In (b) 
the solid, dashed, dotted, and dash-dotted lines correspond to levels 1 - 5 in (a), respectively. 



Z becomes practically impossible. 

3) The principle of compensation of dangerous diagrams was postulated to define the 
coefficients Uj and Vj of the Bogoliubov canonical transformation. This postulation and the 

Q 

variational calculation of dH' /dvj lead to Eq. (19) in [3] for the BCS at T = 0. It is justified 
so long as divergences can be removed from the perturbation expansion of the ground-state 
energy. However, at T 7^ a T-dependent ground state does not exist. Instead, one should 
use the expectation values over the canonical or grand-canonical ensemble . Therefore, 
Eq. (19) of [7] no longer holds at T 7^ since the BCS gap is now defined by Eq. (7) of Q], 
instead of Eq. (3). Fig. IHl clearl3/_^shows how bj 7^ Ej and Cj 7^ at T 7^ 0. This invalidates 
the critics based on Eq. (19) of 7]. 
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In conclusion, the test of 'r\ is inadequate to judge the MBCS apphcabihty because its 
results were obtained either in the T region, where the use of zero-T spectra is no longer valid 
(for ^^°Sn and Ni), or within too limited configuration spaces (the PFM for = ^7 = 10 or 
2 major shells for Ni). Our calculations with a T-dependent spectrum for ^^'^Sn, and within 
extended configuration spaces presented here show that the MBCS is a good approximation 
up to high T even for a system with = 10 particles. 

We thank A. Volya for assistance in the exact solutions of the PFM. 
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